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Abstract 

A  two-dimensional  two-phase  thermal  model  is  presented  for  direct  methanol  fuel  cells  (DMFC),  in  which  the  fuel  and  oxidant  are  fed  in  a  passive 
manner.  The  inherently  coupled  heat  and  mass  transport,  along  with  the  electrochemical  reactions  occurring  in  the  passive  DMFC  is  modeled  based 
on  the  unsaturated  flow  theory  in  porous  media.  The  model  is  solved  numerically  using  a  home- written  computer  code  to  investigate  the  effects  of 
various  operating  and  geometric  design  parameters,  including  methanol  concentration  as  well  as  the  open  ratio  and  channel  and  rib  width  of  the 
current  collectors,  on  cell  performance.  The  numerical  results  show  that  the  cell  performance  increases  with  increasing  methanol  concentration 
from  1.0  to  4.0  M,  due  primarily  to  the  increased  operating  temperature  resulting  from  the  exothermic  reaction  between  the  permeated  methanol 
and  oxygen  on  the  cathode  and  the  increased  mass  transfer  rate  of  methanol.  It  is  also  shown  that  the  cell  performance  upgrades  with  increasing 
the  open  ratio  and  with  decreasing  the  rib  width  as  the  result  of  the  increased  mass  transfer  rate  on  both  the  anode  and  cathode. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  direct  methanol  fuel  cell  (DMFC)  has  been  widely  inves¬ 
tigated  as  an  alternative  power  source  for  powering  various 
potable  electronic  devices  because  they  can  offer  the  promises 
of  higher  energy  efficiency  and  higher  energy  density.  Over  the 
past  decades,  extensive  efforts  have  been  focused  on  the  active- 
feed  DMFC  with  the  fuel  fed  by  a  liquid  pump  and  oxidant 
fed  by  a  gas  compressor  [1-3].  Nevertheless,  these  ancillary 
devices  not  only  make  the  fuel  cell  system  complex  but  also 
lower  the  power  density  because  of  the  parasitic  losses.  For  this 
reason,  the  passive  DMFC  with  neither  liquid  pumps  nor  gas 
compressors  has  been  proposed  and  studied  [4-12].  This  type 
of  passively  operated  fuel  cell  is  appealing  because  it  not  only 
offers  the  advantage  of  a  simpler  and  more  compact  system,  but 
also  eliminates  the  parasitic  power  losses  for  powering  the  ancil¬ 
lary  devices  required  in  the  active  DMFCs.  However,  compared 
to  its  counterpart — the  active  DMFC,  the  performance  of  the 


*  Corresponding  author.  Tel:  +852  2358  8647;  fax:  +852  2358  1543. 
E-mail  address:  metzhao@ust.hk  (T.S.  Zhao). 

0378-7753/$  -  see  front  matter  ©  2007  Elsevier  B.V.  All  rights  reserved, 
doi:  10. 1016/j  .jpowsour. 2007. 09. 086 


passive  DMFC  is  lower  because  of  the  inability  to  handle  the 
excess  liquid  water  on  the  cathode  and  too  much  heat  loss  from 
the  fuel  cell  to  the  ambient,  resulting  in  a  lower  cell  operating 
temperature.  Hence,  the  added  attention  has  to  be  given  to  the 
thermal  and  water  management  for  the  passive  DMFC. 

In  the  DMFC,  the  intrinsically  coupled  physicochemical 
processes  simultaneously  occur,  including  heat  and  mass  trans¬ 
fer,  electrochemical  reactions,  ionic  and  electronic  transfer. 
Since  it  is  difficult  to  experimentally  quantify  the  interrelated 
parameters  governing  the  DMFC,  numerical  modeling  becomes 
essential  for  the  optimization  of  the  cell  design  and  operating 
conditions.  A  number  of  studies  [13-34]  have  been  reported 
on  simulation  of  both  PEMFCs  and  DMFCs.  Most  previous 
models  were  developed  based  on  the  single-phase  flow  the¬ 
ory  [13-26],  which  include  semi-empirical  model  [13]  and 
one-dimensional  and  multi-dimensional  models  [17-25].  How¬ 
ever,  since  the  single-phase  model  cannot  objectively  reflect  the 
two-phase  mass  transport  process  occurring  in  the  DMFC,  as 
revealed  by  the  visualization  study  [3],  the  two-phase  models  for 
the  DMFC  have  recently  been  reported  [27-35].  Murgia  et  al. 
[27]  developed  a  one-dimensional,  two-phase,  multi-component 
steady-state  model  based  on  the  phenomenological  transport 
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Nomenclature 

cj  specific  heat  capacity  of  species  j  (J  mol-1  K-1) 

Cj  concentration  of  species  j  (mol  m-3) 

Dj  diffusivity  of  species  j  (m2  s-1) 

E  thermodynamic  equilibrium  potential  (V) 

F  Faraday’s  constant  (C  mol- 1 ) 

G  Gibbs  free  energy  (J  mol-1) 
h  heat  transfer  coefficient  (W  m-2) 

hw  evaporation  heat  of  water  (J  mol- 1 ) 

H  enthalpy  (J  mol- 1 ) 

i  local  current  density  (A  m-2) 

ip  parasitic  current  density  resulting  from  methanol 
crossover  (A  m-2) 

/ref  exchange  current  density  (A  m-2) 

I  cell  current  density  (Am-2) 

kYg  relative  permeability  of  gas  phase 

kY\  relative  permeability  of  liquid  phase 

kc  condensation  rate  constant  (Is-1) 

kY  evaporation  rate  constant  ( 1  atm- 1  s - 1 ) 

K  permeability  (m2) 

/  channel  and  rib  width  (m) 

Mj  molecular  weight  of  species  j  (kg  mol-1) 
nu2o  electro-osmotic  drag  coefficient  of  water 
n ™  electro-osmotic  drag  coefficient  of  methanol 
Nj  molar  flux  of  species  j  (mol  m-2  s-1) 

P  atmosphere  pressure  (Pa) 

Pc  pressure  (Pa) 

Pg  gas  pressure  (Pa) 

P\  liquid  pressure  (Pa) 

q  heat  generation  rate  (W  m-2) 

R  interfacial  transfer  rate  or  gas  constant 

(molm-3  s-1  or  J  mol-1  K-1) 

^contact  contact  resistance  (£2  m2) 
s  liquid  saturation  in  porous  media 

T  temperature  (K) 

u  superficial  velocity  vector  (ms-1) 

V  potential  (V) 

v  coordinate 

y  coordinate  or  mole  fraction 

Y*  dimensionless  coordinate 

Greek  letters 

aa  charge  transfer  coefficient  at  anode 

ac  charge  transfer  coefficient  at  cathode 
8  thickness  of  the  membrane  (m) 

y)  overpotential  (V) 

0C  contact  angle  (°) 

k  ionic  conductivity  of  membrane  (1  m-1) 

/x  viscosity  (kgm-1  s-1) 

p  density  (kgm-3) 

a  surface  tension  (N  m- 1 ) 


Superscripts  and  subscripts 
a  anode 

acl  anode  catalyst  layer 

agdl  anode  gas  diffusion  layer 

c  cathode 

cell  cell 

ccl  cathode  catalyst  layer 

cgdl  cathode  gas  diffusion  layer 

ccc  cathode  current  collector 

eff  effect  value 

g  gas  phase 

H2O  water 

1  liquid  phase 

m  methanol 

mem  polymer  electrolyte  membrane 

O2  oxygen 

ref  reference  value 

sat  saturated 

v  vapor 


equations.  This  model  introduced  a  Gaussian  function  to  approx¬ 
imately  account  for  the  influence  of  the  capillary  pressure  on 
the  effective  gas  porosity.  Wang  and  Wang  [28]  developed  a 
model  for  the  DMFC  based  on  the  mixture  multi-phase  flow 
model.  Divisek  et  al.  [30]  presented  a  two-dimensional,  two- 
phase  model  for  the  DMFC  with  the  computational  domain  of  a 
cross-section  perpendicular  to  the  flow  direction  at  a  certain  posi¬ 
tion,  and  the  model  was  restricted  to  the  porous  electrodes  and 
membrane  of  the  cell.  In  their  work,  the  presence  of  hydrophilic 
and  hydrophobic  pores  was  taken  into  account  by  introducing  a 
supposed  relationship  between  capillary  pressure  and  saturation 
which  is  different  from  the  Levereete  function  as  used  in  [28]. 
In  addition,  the  two-phase  mass  transport  phenomena  in  PEM- 
FCs  have  been  modeled  based  on  the  unsaturated  flow  theory  in 
porous  media  (UFT)  [31-34]. 

In  addition  to  modeling  of  active  DMFCs  and  PEMFCs  men¬ 
tioned  above,  the  mathematical  models  for  the  air-breathing 
PEMFCs  and  passive  DMFCs  have  also  been  reported  [36-48]. 
Ziegler  et  al.  [42]  developed  a  single-phase  two-dimensional 
non-isothermal  model  for  an  air-breathing  PEMFC  and  inves¬ 
tigated  the  effect  of  the  geometric  design  of  the  cathode  end 
plate  on  the  gas,  potential  and  cell  temperature  distribution. 
Hwang  [44]  presented  a  three-dimensional  model  to  study 
the  species  transport  and  electrochemical  characteristics  in 
a  free-breathing  cathode  of  a  planar  fuel  cell,  in  which  the 
Brinkman  extension  to  the  Darcy’s  law  was  employed  to  describe 
the  gas-mixture  flow  in  the  porous  cathode.  The  effects  of 
breathing-hole  arrangements  and  open-area  ratio  of  the  cur¬ 
rent  collector  on  the  electrochemical  performance  were  also 
examined.  Recently,  Rice  and  Faghri  [48]  developed  a  tran¬ 
sient  two-dimensional,  multi-phase,  multi-component  model  for 
the  DMFC  with  a  passive  fuel  delivery  system.  The  passive 
delivery  system  utilized  a  porous  medium  to  passively  deliver 
methanol  to  the  fuel  cell  while  controlling  the  concentration 
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of  methanol  at  the  anode  side  to  limit  the  amount  of  methanol 
crossover. 

Our  literature  review  indicates  that  the  two-dimensional  two- 
phase  model  with  the  heat  transfer  effect  for  the  passive  DMFC 
has  not  been  reported  yet.  However,  as  mentioned  above,  the 
thermal  and  water  management  is  crucial  for  the  passive  DMFC. 
Hence,  the  two-phase  transport  in  the  passive  DMFC  with  taking 
account  of  the  heat  transfer  effect  is  essential  for  better  under¬ 
standing  of  transport  phenomena  in  this  type  of  fuel  cell.  In 
this  work,  we  developed  a  two-dimensional  two-phase  thermal 
model  for  the  passive  DMFC.  The  two-phase  transport  in  the 
anode  and  cathode  porous  region  (gas  diffusion  layer)  is  formu¬ 
lated  based  on  the  unsaturated  flow  theory  in  porous  media.  With 
this  model,  the  effects  of  species  distribution,  current  density  dis¬ 
tribution  and  temperature  distribution  as  well  as  the  effects  of 
the  various  operating  and  geometric  design  parameters  on  the 
cell  performance  are  studied. 

2.  Mathematical  model 

Consider  the  passive  DMFC  sketched  in  Fig.  la,  which  con¬ 
sists  a  fuel  tank,  an  anode  current  collector  (ACC),  an  anode 
gas  diffusion  layer  (AGDL),  an  anode  catalyst  layer  (ACL),  a 
polymer  electrolyte  membrane  (MEM),  a  cathode  catalyst  layer 
(CCL),  a  cathode  gas  diffusion  layer  (CGDL),  and  a  cathode  cur¬ 
rent  collector  (CCC).  On  the  anode,  methanol  is  transported  from 
the  fuel  tank  to  the  ACL,  where  part  of  methanol  is  electrochem- 
ically  oxidized  to  form  CO2,  while  the  remainder  permeates  the 
membrane  and  reacts  with  oxygen  on  the  cathode.  The  produced 
gas  CO2  transfers  back  to  the  fuel  tank  and  is  emitted  to  the  ambi- 


Fig.  L  (a)  Schematic  of  the  passive  DMFC;  (b)  the  computation  domain. 


ent.  On  the  cathode,  oxygen  is  passively  taken  from  the  ambient 
air  to  the  CCL,  where  it  is  reduced  to  form  water  and  heat.  The 
generated  water  moves  out  through  the  CGDL  by  the  capillary 
force.  The  heat  generated  in  the  fuel  cell  is  lost  from  the  CCC 
to  the  ambient  by  natural  convection. 

To  save  the  computational  expense,  a  unit  region,  as  illus¬ 
trated  in  Fig.  lb,  can  be  taken  from  the  entire  cell  as  the 
computational  domain,  in  which  the  components  in  the  hori¬ 
zontal  direction  remain  the  same  as  in  Fig.  la,  but  in  the  vertical 
direction  only  the  half  channel  and  half  rib  width  are  included 
because  of  the  lateral  symmetry  with  respect  to  the  middle  point 
of  channel  width.  Note  that  the  v-axis  and  y-axis  origins  are  set 
at  the  outer  surface  of  the  AGDL  and  at  the  middle  of  channel 
width,  respectively.  For  convenience  of  description,  each  bound¬ 
ary  and  interface  are  designated  with  Roman  numerals  from  I 
to  X.  The  model  is  formulated  based  on  the  following  general 
simplifications  and  assumptions: 

(1)  The  fuel  cell  is  assumed  to  operate  under  steady-state  con¬ 
ditions. 

(2)  The  methanol  concentration  in  the  fuel  tank  remains  con¬ 
stant  during  the  fuel  cell  discharging. 

(3)  The  fuel  tank,  ACC,  and  AGDL  are  well  insulated  from 
the  ambient.  No  heat  is  lost  from  these  components  to  the 
ambient.  Thus,  the  temperature  at  each  of  these  components 
is  the  same  as  that  at  ACL. 

(4)  Considering  that  ACL  and  CCL  are  much  thinner  than 
AGDL,  CGDL  and  MEM,  they  are  treated  as  an  interface 
and  the  temperature  and  concentration  distributions  are  uni¬ 
form. 

(5)  On  the  cathode,  the  heat  is  transferred  from  the  CCC  to  the 
ambient  by  natural  convection. 

(6)  Since  the  reaction  rate  of  methanol  on  the  cathode  is  rather 
fast,  the  methanol  in  the  CCL  is  assumed  to  be  completely 
depleted. 

(7)  Since  the  variation  in  the  fluid  viscosity  and  density  with 
temperature  is  relatively  small  in  the  calculated  range,  these 
properties  are  considered  to  be  constant. 

(8)  Considering  only  CO2  in  the  anode  gas  phase,  the  methanol 
vapor  and  water  vapor  transport  is  ignored. 

2.7.  Two-phase  mass  transport  in  the  AGDL 

We  first  consider  the  two-phase  mass  transport  in  the  AGDL. 
Based  on  unsaturated  flow  theory  in  porous  media,  on  the  anode, 
the  liquid  pressure  is  assumed  to  be  constant.  Thus,  the  transfer 
of  methanol  from  the  fuel  tank  to  the  ACL  relies  on  diffusion 
only,  where  part  of  methanol  is  electrochemically  oxidized  to 
produce  gas  CO2,  electrons  and  protons,  while  the  remainder 
permeates  the  membrane  and  arrives  in  the  CCL.  The  produced 
gas  CO2  transfers  back  to  the  fuel  tank  by  the  capillary  force. 
This  two-phase  counter-flow  mass  transport  in  the  AGDL  can 
be  modeled  by: 

V  •  (£>^VC^gdl)  =  0  (liquid phase), 


(1) 
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and 


(gas  phase), 


(2) 


where  D^,  pg  and  Mg  represent  the  effective  diffusivity  of 
methanol  in  the  AGDL,  the  density  and  molar  weight  of  gas 
CO2,  respectively.  In  Eq.  (2),  ug  is  the  superficial  velocity  vec¬ 
tor  based  on  the  total  cross-sectional  area  of  fluids  and  porous 
medium,  which  can  be  obtained  from  Darcy’s  law: 


u 


g  ~ 


Kkvg 

H 


v^g 


(3) 


with  K ,  kYg  and  /xg  denoting  the  absolute  permeability  of  the 
AGDL,  the  relative  permeability  of  the  gas  phase  and  the  vis¬ 
cosity  of  the  gas  phase,  respectively.  The  pressure  difference 
between  the  gas  phase  and  liquid  phase  is  related  to  the  capillary 
pressure: 


(  s  \  °-5 

Pc  =  Pg  ~  P\  =  cr  cos(Pc)  ^  —  j  J(s ), 


(4) 


where  J(s)  represents  the  Levertte  function  with  s  being  as  liq¬ 
uid  saturation  in  a  porous  medium  (see  Table  1).  Based  on  the 
assumption  of  the  constant  liquid  pressure  in  the  AGDL,  we 
combine  Eqs.  (2)-(4)  to  obtain: 


V  • 


f  Pg  Kkrg  ^ 
V  Mg  Pg 


0.5  \ 

V/(j)  J 


=  0. 


(5) 


2.2.  Two -phase  mass  transport  in  the  CGDL 


gas  phase  is  constant,  we  can  treat  this  mass  transport  process 
in  the  gas  phase  to  be  a  diffusion  problem.  As  a  result,  this  two- 
phase  counter-flow  mass  transport  in  the  CGDL  can  be  modeled 
with  the  following  equations: 

-v  •  (-D^vcgf )  =  0  (gas  phase),  (6) 

-V  •  (-Df' VCcfdl)  -  =  0  (gas  phase),  (7) 

and 

— V  •  +  Rw  =  0  (liquid phase),  (8) 

where  Dq^  and  stand  for  the  effective  diffusivity  of  the 
oxygen  and  water  vapor;  p\  and  M\  denote  the  density  and  molar 
weight  of  the  water;  Rw  is  the  interfacial  transfer  rate  of  water 
between  liquid  and  vapor  water  (see  Table  1);  ui  is  the  superficial 
velocity  vector  based  on  the  total  cross-sectional  area  of  fluids 
and  porous  medium,  which  can  be  obtained  from  Darcy’s  law: 


ui  =  - 


Kk d 
P\ 


V/?i, 


(9) 


where  kr\  is  the  relative  permeability  of  the  liquid  phase,  p\  is 
the  viscosity  of  the  liquid  water,  p\  is  the  pressure  of  the  liquid 
phase.  Similarly,  according  to  the  definition  of  the  capillary  force 
shown  in  Eq.  (4)  and  the  assumption  of  the  constant  gas  pressure 
on  the  cathode,  we  can  obtain: 


V-  - 


Pi  Kkr\ 
M\  /xi 


crcos(0c)(-^)  V ./(.?) 


+  Rw  —  0. 


(10) 


We  now  turn  our  attention  to  the  two-phase  mass  transport  in 
the  cathode  gas  diffusion  layer  (CGDL).  Oxygen  transfers  to  the 
cathode  catalyst  layer  (CCL),  where  part  of  the  oxygen  is  elec- 
trochemically  reduced  to  form  liquid  water,  while  the  remainder 
directly  reacts  with  the  permeated  methanol  to  produce  heat  and 
liquid  water  and  CO2.  The  liquid  water  generated  in  the  CCL 
and  that  permeated  from  the  anode  move  to  the  ambient  by  the 
capillary  force.  It  should  be  noted  that  this  process  is  coupled 
with  the  liquid  water  evaporation.  Assuming  that  the  pressure  of 


2.3.  Mass  transport  in  the  membrane 


We  now  consider  the  transport  of  methanol  and  water  through 
the  membrane.  Methanol  permeates  the  membrane  as  the  result 
of  the  concentration  gradient  and  electro-osmotic  drag  by  the 
proton  transfer.  Accordingly,  the  flux  of  methanol  through  the 
membrane,  Across,  can  be  determined  from: 


N cross  —  T) 


,eff 


dC 


m,mem  .  m  l 

—  +"‘T 


(ii) 


Table  1 

Constitutive  relations  and  definitions 


Parameters 

Expressions 

Capillary 

Pc  =  Pg 

—  p\  =  0  cos  6c(£/K)05  J(s) 

pressure 

m  =  < 

f  1.417(1  -  s)  -  2.120(1  -  s)2  +  1.263(1  -  s)3  0  <  0C  <  90° 

^  1.417^  -  2.12052  +  1.2635s  90°  <  <9C  <  180° 

Relative 

kr\  =  S3 

liquid 

permeability 

pg  =  (1 

-  s)3  gas 

Effective 

Df  =  Di£15(l  -  s)1'5  i:02,  vapor  CGDL 

diffusion 

Df  =  Die'-V'5  i  :  methanol  AGDL 

coefficients  of 

species 

The  saturation  pressure  of  vapor 

-2.1794  +  0.02953(7  -  273) 

logio  pf  =  -9.1837  x  10“5(T  -  273)2  atm 

( 

+1.4454  x  107(7  -  273)3 

ke  Mh  0  (')’vPs  P ^  ywPs  <  ^ 

ky  RT  (yvPg  pf)  y^Pi  >  pf 

Interfacial  transfer  rate  of  water  between  liquid  and  vapor  [34] 

Rw=  < 
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where  D^fmem  is  the  effective  diffusivity  of  methanol  in  the 
membrane,  F  is  the  Faraday’s  constant  (96,485  C  mol-1),  n ™ 
is  the  electro-osmotic  drag  coefficient  of  methanol.  The  trans¬ 
port  of  liquid  water  through  the  membrane  is  also  driven  by  the 
electro-osmotic  drag  and  concentration  gradient.  However,  since 
the  concentration  difference  of  water  across  the  membrane  is 
rather  small,  the  flux  of  water  can  be  determined  by  considering 
the  electro-osmotic  drag  only,  i.e.: 

i 

A^h2o=^h2o—  (12) 

r 

with  ftH2o  denoting  electro-osmotic  drag  coefficient  of  water. 

2.4.  Heat  transfer 


We  now  consider  the  heat  transfer  process  in  the  passive 
DMFC.  On  the  anode,  the  heat  generated  by  the  electrochemical 
reaction  in  the  ACL,  is  given  by: 


AH,-  AGa\ 
6F  )  ’ 


(13) 


where  the  first  term  represents  the  heat  generation  due  to  the  acti¬ 
vation  and  mass  transfer  overpotentials  on  the  anode,  the  second 
term  represents  the  entropy  change  of  the  anodic  electrochem¬ 
ical  reaction,  with  AH,  denoting  the  anodic  reaction  enthalpy 
and  AGa  the  Gibbs  free  energy. 

The  heat  generation  in  the  CCL,  qcc i,  can  be  determined  from: 


qcd  =  (i  +  ip)0c  ~  (i  +  *p)- 


AHC  -  AGC 


AH,  -  AG, 


4  F 


l  D 


6F 


(14) 


we  assume  the  MOR  is  a  first-order  reaction.  Thus,  the  current 
density  can  be  determined  by  the  simplified  Tafel  equation: 


•ref  ^m,acl^ 


/^ref 


exp 


(16) 


where  and  Cjjf  represent,  respectively,  the  reference 
exchange  current  density  on  the  anode  and  the  reference  con¬ 
centration  of  methanol,  Cm;acl  and  s  stand  for  the  methanol 
concentration  and  the  liquid  saturation  at  the  ACL. 

On  the  cathode,  to  account  for  the  effect  of  methanol 
crossover  on  the  cathode  overpotential,  we  assume  that  the 
methanol  permeated  from  the  anode  completely  and  electro- 
chemically  reacts  with  the  oxygen  on  the  cathode.  Based  on  this 
assumption,  the  internal  current,  /p,  due  to  the  methanol  oxida¬ 
tion  on  the  cathode  can  be  obtained  based  on  the  permeation  flux 
of  methanol  [28]: 


-P 


(17) 


It  follows  from  Eqs.  (16)  and  (17)  that  the  cathode  overpo¬ 
tential  taking  account  of  the  mixed  potential  can  be  determined 
from: 


•ref  6^02,ccl(l  s)  (  \ 

r= cjg  “p  (  MS*,)  '  <18) 

where  {q  and  Cq1  represent,  respectively,  the  reference 
exchange  current  density  on  the  cathode  and  the  reference 
concentration  of  oxygen,  Co2,ccl  and  s  denote  the  oxygen  con¬ 
centration  and  liquid  saturation  at  the  CCL. 


where  the  first  term  represents  the  heat  generation  due  to  activa¬ 
tion  and  mass  transfer  overpotentials  and  mixed  potential  caused 
by  methanol  crossover  on  the  cathode,  the  second  term  accounts 
for  the  entropic  loss,  with  A Hc  denoting  the  cathodic  reaction 
enthalpy  and  A  Gc  the  Gibbs  free  energy,  and  the  third  term 
denotes  the  entropy  change  of  methanol  oxidation  reaction  on 
the  cathode  due  to  methanol  crossover. 

Since  the  fuel  tank,  ACC,  and  AGDL  are  assumed  to  be  well 
insulated  from  the  ambient,  the  heat  generated  within  the  cell  is 
lost  to  the  ambient  from  the  cathode.  The  heat  transfer  process 
involves  three  components:  the  membrane,  CGDL  and  CCC. 
Since  the  flow  velocity  in  the  passive  DMFC  is  rather  small,  the 
heat  transfer  process  through  the  membrane,  CGDL  and  CCC 
can  be  modeled  as  a  heat-conduction  problem  as  follows: 

V  •  (-A ,*vr)  =  0,  i  :  MEM,  CCC,  (15a) 

and 

V  •  (-A ,*vr)  =  hvRw,  i  :  CGDL,  (15b) 

where  A ,*■  is  the  effective  thermal  conductivity  of  each  compo¬ 
nent,  hy  is  the  latent  heat  of  vaporization  of  liquid  water,  Rw  is 
the  water  evaporation  rate. 

2.5.  Electrochemical  kinetics 

On  the  anode,  the  methanol  oxidation  reaction  is  a  compli¬ 
cated  multi-step  electrochemical  reaction.  For  simplification, 


2.6.  Boundary  conditions 

At  the  outer  surface  of  the  AGDL  and  CGDL,  the  methanol 
concentration,  oxygen  concentration  and  water  vapor  concentra¬ 
tion  are  given  based  on  the  reactant  supplying  conditions.  The 
internal  interfacial  conditions  are  given  based  on  the  continuity 
and  mass/species  flux  balance  to  satisfy  the  general  mass  and 
species  conservation  of  the  entire  cell.  The  boundary  conditions 
at  each  interface  shown  in  Fig.  lb  are  detailed  in  Table  2. 

2. 7.  Cell  performance 


The  protons  and  electrons  produced  from  the  MOR  in  the 
ACL  transfer  to  the  cathode  through  the  membrane  and  the  exter¬ 
nal  circuit,  respectively.  The  mean  cell  current  density  can  be 
calculated  by: 


j  pyy 

i 


(19) 


Similarly,  on  the  cathode,  the  mean  parasitic  current  density 
generated  by  the  cathode  methanol  oxidation  can  be  determined 
from: 


h  = 


/  h  d>: 

z 


(20) 


With  the  methanol/oxygen  concentration  at  the  catalyst  layers 
and  the  temperature  distributions,  and  the  anodic  and  cathodic 
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Table  2 

Boundary  conditions 

Boundary  conditions 


I 

II 

III 

IV 

V 

VI 

VII 

VIII 

IX 
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dx  dx 
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6  F 

AH* -AG*' 


Kpgkrg  V7  _  i  _ 
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Mgfig  6  F 


‘l* — 

i  3 
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i  i 


Nu2o  ■ 


2NCTOSs  +  nn2o  ^  T 


AHC-AGC 


IF' 


Co2  —  C02,00?  Cv  —  Cv,< 

/z(T  -  Too) 

&Cp2  _  Q  dCy  _  Q 
dx  dx 

d  T 

k-ccc  —  h(T  —  Tqq) 

dy 
d  T 

—7Ccc  —  h(T  —  Too) 

dx 

d^m  _  Q  dCp2  _  ^ 
dy  ’  dy 

d^m  _  q  dCp2  _  g 

dy  ’  dy 


qcc\  =  (i  +  ip)rjc  ~  i~ 

nF 
d  T 

x),  s  =  0.05,  ^cgdl — ;  — 


dx 


ds  dT  _  dr 

—  —  0,  rcgdi  —  —  rccc  — 

dx  dx  dx 


dCy  ds  dr 

-7^=0,  —  =  0,  —  =  0 

dy  dy  dy 

dCv  _  ds  dr 

-7^=0,  —  =  0,  —  =  0 

dy  dy  dy 


overpotentials  obtained  from  the  equations  presented  above,  we 
can  assess  the  cell  performance  based  on: 

5 

^cell  =  ^cell  —  Oa.  ~~  0c  ~  I  ^contact?  (21) 

K 

where  £cell  is  the  thermodynamic  equilibrium  potentials  of  the 
fuel  cell  and  is  a  function  of  temperature  and  pressure,  Contact 
is  the  contact  resistance  of  the  fuel  cell.  The  thermodynamic 
equilibrium  potentials  of  the  fuel  cell  can  be  calculated  by: 

Ecell  =  £c%  +  A  T  (jEj  ,  (22) 

where  is  the  open  circuit  voltage  at  T=298K  and 

(dE/dT)  represents  the  change  rate  of  electromotive  force  of 
-1.4  x  10-4  VK-1  [49]. 

3.  Results  and  discussion 

Using  a  self- written  code  based  on  the  SIMPLE  algorithm,  all 
the  variables  are  iteratively  solved  based  on  the  physicochemical 
properties  given  in  Table  3  until  the  pre-set  convergent  criteri- 
ons  are  met.  With  this  model,  the  effects  of  several  parameters 
including  the  methanol  concentration  and  collector  design  on 
cell  performance  were  investigated.  Details  are  presented  below. 

3.1.  Effect  of  methanol  concentration 

Fig.  2  shows  the  polarization  curves  when  the  fuel  cell  oper¬ 
ates  with  1-4  M  methanol  solutions,  respectively.  It  can  be  seen 
from  this  figure  that  the  cell  performance  upgrades  progressively 
with  an  increase  in  methanol  concentration.  The  increased  cell 
performance  as  the  result  of  increasing  methanol  concentration 
can  be  explained  as  follows.  Firstly,  the  mass  transfer  rate  of 


DMFC. 

methanol  at  the  anode  increases  with  methanol  concentration, 
reducing  the  voltage  loss  on  the  anode  due  to  the  mass  transfer 
polarization  and  thereby  tending  to  yield  better  cell  performance. 
Secondly,  Fig.  3  also  shows  that  the  rate  of  methanol  crossover 
increases  with  methanol  concentration.  The  increased  rate  of 
methanol  crossover  results  in  a  larger  mixed  potential,  which 
tends  to  reduce  the  cell  voltage.  Meanwhile,  it  is  also  seen  that 
the  rate  of  methanol  crossover  decreases  with  the  current  den¬ 
sity  as  a  result  of  the  reduced  methanol  concentration  at  the 
anode  electrode/membrane  interface.  Thirdly,  the  increased  rate 
of  methanol  crossover  as  the  result  of  higher  methanol  con¬ 
centration  will  also  lead  to  a  higher  rate  of  methanol  oxidation 
on  the  cathode,  which  in  turn  results  in  a  higher  heat  gener¬ 
ation  rate.  As  a  result,  the  cell  operating  temperature  virtually 
increases  with  methanol  concentration,  as  evidenced  from  Fig.  4, 
which  shows  the  variation  in  the  temperature  at  the  ACL  with 
current  density  for  different  methanol  concentrations.  In  partic¬ 
ular,  it  is  seen  that  the  temperature  difference  between  different 
methanol-concentration  operations  at  low  current  densities  and 
open  circuit  condition  (OCC)  is  much  larger  than  that  at  mid¬ 
dle  and  high  current  densities.  This  is  because  the  methanol 
crossover  is  more  serious  at  low  current  densities  and  OCC, 
leading  to  a  higher  heat  generation  rate  and  thereby  larger  tem¬ 
perature  difference.  The  increased  operating  temperature  fasters 


Fig.  3.  Variations  in  the  methanol  crossover  current  with  current  density  for 
different  methanol  concentrations. 
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Table  3 

Physicochemical  properties 

Parameter,  symbol  (unit)  Value  Refs. 


Liquid  methanol  enthalpy  of  formation,  A Hm  (J  mol-1) 

Liquid  water  enthalpy  of  formation,  AHW  (J  mol-1) 

Carbon  dioxide  enthalpy  of  formation,  AHqo2  (J  mol-1) 

Liquid  methanol  Gibbs  free  energy,  AGm  (J  mol-1) 

Liquid  water  Gibbs  free  energy,  A Gw  (J  mol-1) 

Carbon  dioxide  Gibbs  free  energy,  AGco2  (J  mol-1) 

Liquid  methanol  specific  heat  capacity,  cm  (J  mol-1  K-1) 

Liquid  water  specific  heat  capacity,  cw  (J  mol-1  K-1) 

Carbon  dioxide  specific  heat  capacity,  cqo2  (J  mol-1  K-1) 
Oxygen  specific  heat  capacity,  Co2  (Jmol-1  K-1) 

Evaporation  heat  of  water,  hw  (J  mol-1) 

Pressure  of  air  in  cathode,  Pg  (Pa) 

Thickness  of  ACC  and  CCC,  /ccc  (m) 

Thickness  of  AGDL  and  CGDL,  /gdi  (m) 

Electro-osmotic  drag  coefficient  of  water,  nu2o 
Electro-osmotic  drag  coefficient  of  methanol,  nJJ1 
Proton  conductivity  in  membrane,  k  (1  £2-1  m-1) 

Porosity  of  AGDL  and  CGDL,  s 
Permeability  in  AGDL  and  CGDL,  K  (m2) 

Thermal  conductivity  of  membrane,  A.mem  (W  m-1  K-1) 

Thermal  conductivity  of  CGDL,  kcgdi  (W  m-1  K-1) 

Thermal  conductivity  of  CCC,  7CCC  (W  m-1  K-1) 

Heat  transfer  coefficient,  h  (W  m-1  K-1) 

Diffusion  coefficient  of  methanol  in  membrane,  D^JTmem  (m2  s-1) 
Diffusivity  of  methanol,  Dm  (m2  s-1) 

Diffusivity  of  oxygen,  Do2  (m2  s-1) 

Diffusivity  of  water  vapor,  Dv  (m2  s-1) 

Density  of  gas  phase,  pg  (kgm-3) 

Density  of  liquid  phase,  p\  (kg  m-3) 

Viscosity  of  gas  phase,  /zg  (kgm-1  s-1) 

Viscosity  of  liquid  phase,  (kgm-1  s-1) 

Surface  tension,  o  (N  m-1) 

Transfer  coefficient  of  anode,  aa 
Transfer  coefficient  of  anode,  ac 
Reference  concentration  of  methanol,  Cjjf  (molm-3) 

Reference  concentration  of  oxygen,  Cq1  (molm-3) 

Evaporation  rate  constant,  ke  (1  atm-1  s-1) 

Condensation  rate  constant,  kc  (Is-1) 

Reference  exchange  current  density  on  anode,  iJJf  (Am-2) 
Reference  exchange  current  density  on  cathode,  z'q1  (Am-2) 
Contact  resistance,  ^contact  cm2) 


-238.66  x  103 
-285.83  x  103 
-393.51  x  103 
-166.27  x  103 
-237.08  x  103 
-394  x  103 
80.96 
75.24 

36.9 
39.44 

44.86  x  103 
1  x  105 
0.0005 


0.0003  [50] 

2.9  exp(  1029/(1/3 3 3  -  1/7))  [23] 

%  x  nu2o  [28] 

7.3  exp(1268(l/298  -  1/7))  [14] 

0.6  [33] 

1.2  xlO-12  [34] 

0.21  [24] 

1.6  [24] 

16  [24] 

20  [42] 

4.9  x  10-10  exp(2436(l/333  -  1/7))  [14] 

2.8  x  10-9  exp  (2436(1/333  -  1/7))  [14] 

0.1775  x  10-4  x  (77273.15)1-823  [33] 

0.256  x  10-4  x  (77273. 15)2'334  [33] 


1.98 
1  x  103 
14.96  x  10-6 
0.9  x  10-3 
0.0625 


0.35  [24] 

0.8  [24] 

4  x  103  [24] 

0.21  x  Pg/RT 

1  [33] 

100  [33] 

94.25  exp(35570/R*(l/353  -  1/7))  [28] 

0.04222  exp(73200/R*(l/353  -  1/7))  [28] 

0.8  [23] 


Fig.  4.  Variations  in  the  operating  temperature  at  the  ACL  with  current  density 
for  different  methanol  concentrations. 


the  electrochemical  kinetics  on  both  the  anode  and  cathode,  pro¬ 
ducing  a  higher  voltage  at  a  given  current  density.  Therefore,  the 
improved  performance  with  increasing  methanol  concentration 
shown  in  Fig.  2  can  be  attributed  to  the  increased  cell  operating 
temperature  and  the  increased  mass  transfer  rate  of  methanol, 
although  the  mixed  potential  on  the  cathode  also  increases  with 
methanol  concentration.  The  variation  in  the  cell  performance 
with  methanol  concentration  is  consistent  with  that  reported 
elsewhere  [24] .  It  is  also  noted  from  Fig.  4  that  for  each  methanol 
concentration,  the  temperature  at  the  ACL  increases  with  current 
density.  This  is  because  the  heat  generation  rate  due  to  the  kinetic 
and  mass  transfer  polarization  increases  with  current  density. 

3.2.  Effect  of  current-collector  open  ratio 

The  effect  of  the  current-collector  open  ratio  was  studied  by 
equipping  the  DMFC  with  identical  current  collectors  at  both 
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the  anode  and  cathode.  The  open  ratio  was  changed  by  fixing 
the  sum  of  channel  and  rib  width  at  2.0  mm.  Fig.  5  shows  that 
the  cell  performance  increases  with  the  open  ratio,  particularly 
at  high  current  densities.  The  improved  cell  performance  as  the 
result  of  increasing  the  open  ratio  can  be  primarily  attributed 
to  the  enhanced  mass  transfer  rates  of  both  methanol  and  oxy¬ 
gen.  The  reasons  leading  to  the  increased  mass  transfer  rates  are 
as  follows.  Firstly,  the  specific  area  of  mass  transfer  increases 
with  the  open  ratio,  yielding  higher  mass  transfer  rates  of  both 
methanol  and  oxygen.  Secondly,  the  rib  width  decreases  with  the 
open  ratio,  providing  a  shorter  distance  of  mass  transfer  from  the 


Fig.  6.  Variations  in  the  operating  temperature  at  the  ACL  with  current  density 
for  different  open  ratios. 

channel  region  to  the  rib  region  in  the  GDL  and  thereby  result¬ 
ing  in  the  higher  methanol  and  oxygen  concentration  under  the 
ribs.  Thirdly,  the  gas  CO2  on  the  anode  and  liquid  water  on  the 
cathode  are  more  easily  to  be  removed  from  the  rib  region  to  the 
outside  as  the  open  ratio  increases.  Subsequently,  more  active 
sites  can  be  used  for  the  electrochemical  reaction  in  the  catalyst 
layer,  leading  to  the  improved  cell  performance. 

The  variations  in  the  cell  temperature  at  the  ACL  with  the 
current  density  for  different  open  ratios  are  shown  in  Fig.  6. 
For  each  open  ratio,  the  cell  temperature  increases  with  the  cur- 


Fig.  7.  Distribution  of  (a)  methanol  concentration  at  the  ACL,  (b)  current  density,  (c)  gas  saturation  at  the  ACL,  (d)  liquid  saturation  at  the  CCL  with  different  open 
ratios. 
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rent  density,  as  the  heat  generation  rate  due  to  the  kinetic  and 
mass  transfer  polarization  increases  with  the  current  density.  It 
is  interesting  to  notice  that  at  low  current  densities,  the  larger 
open  ratio  yields  a  slightly  higher  cell  operating  temperature 
than  does  the  smaller  one,  but  at  high  current  densities  the  trend 
reverses.  This  is  because  the  larger  open  ratio  yields  a  higher 
rate  of  methanol  crossover  at  low  current  densities,  resulting  in 
a  higher  mixed  potential  on  the  cathode,  which  in  turn  causes  a 
higher  heat  generation  rate  on  the  cathode.  As  a  result,  the  larger 
open  ratio  yields  a  slightly  higher  cell  temperature  at  low  cur¬ 
rent  densities.  As  the  current  density  increases,  the  overpotential 
due  to  the  mass  transfer  becomes  more  dominant.  The  increased 
mass  transfer  loss  with  the  smaller  open  ratio  produces  more 
heat  at  high  current  densities,  leading  to  a  higher  cell  tempera¬ 
ture.  Moreover,  the  increased  rib  region  as  the  result  of  reduction 
in  the  open  ratio  also  increases  the  heat  transfer  resistance  from 
the  cell  to  the  ambient,  reducing  the  heat  loss.  This  also  con¬ 
tributes  to  the  higher  cell  temperature  for  the  smaller  open  ratio 
at  high  current  densities.  It  should  be  pointed  that  since  the  cell 
temperature  change  with  the  open  ratio  is  rather  small,  it  has  lit¬ 
tle  effect  on  the  cell  performance.  Therefore,  the  improved  cell 
performance  as  the  result  of  increasing  the  open  ratio  shown  in 
Fig.  5  can  be  mainly  attributed  to  the  increased  rates  of  methanol 
transfer  on  the  anode  and  oxygen  transfer  on  the  cathode. 

Fig.  7  shows  the  distributions  of  methanol  concentration  at 
the  ACL,  current  density,  gas  saturation  at  the  ACL,  and  liquid 
saturation  at  the  CCL  for  different  open  ratios  at  a  given  anode 
overpotential  of  0.5  V.  It  is  seen  from  Fig.  7a  that  a  more  uniform 
distribution  of  methanol  concentration  is  achieved  with  increas¬ 
ing  the  open  ratio  as  the  result  of  the  enhanced  mass  transfer 
mentioned  earlier,  which  in  turn  leads  to  higher  and  more  uni¬ 
form  distribution  of  current  density,  as  shown  in  Fig.  7b.  These 
results  further  confirm  that  the  improved  cell  performance  shown 
in  Fig.  5  is  primarily  due  to  the  enhanced  mass  transfer  rates 
resulting  from  the  large  open  ratio.  Meanwhile,  a  higher  current 
density  produces  more  gas  CO2  on  the  anode  and  liquid  water  on 
the  cathode.  As  a  result,  the  larger  open  ratio  yields  higher  gas 
saturation  at  the  ACL  and  liquid  saturation  at  CCL,  as  shown  in 
Fig.  7c  and  d.  It  is  also  found  that  since  the  gas  CO2  is  difficult 
to  transport  from  the  rib  region  to  the  outside  due  to  the  larger 
mass  transfer  resistance,  leading  to  higher  gas  saturation  under 
the  rib.  However,  this  is  not  the  case  for  the  liquid  saturation  at 
the  CCL.  The  liquid  saturation  under  the  rib  is  lower  than  that 
under  the  channel  region  because  of  the  lower  water  generation 
rate. 

The  cell  temperature  distributions  in  the  x  and  y  directions 
for  different  open  ratios  at  a  given  anode  overpotential  of  0.5  V 
are  presented  in  Fig.  8.  It  is  seen  from  this  figure  that  a  larger 
open  ratio  exhibits  a  higher  cell  temperature  because  a  higher 
current  density  leads  to  a  higher  heat  generation  rate.  It  is  inter¬ 
esting  to  notice  that  the  variation  trend  in  the  cell  temperature 
shown  in  Fig.  8a  is  similar  to  the  current  density  distribution 
shown  in  Fig.  7b.  This  is  because  the  heat  generation  rate  is 
proportional  to  the  current  density.  The  higher  current  density 
at  the  channel  region  results  in  the  higher  heat  generation  rate 
and  thus  higher  cell  temperature.  Since  the  heat  is  generated  at 
the  ACL  and  CCL  and  released  from  the  cathode  to  the  ambient, 


Fig.  8.  Distribution  of  (a)  cell  temperature  in  y  direction,  (b)  cell  temperature  in 
x  direction  with  different  open  ratios. 

the  cell  temperature  decreases  along  the  x  direction,  as  shown  in 
Fig.  8b.  Particularly,  the  larger  open  ratio  shows  a  little  higher 
temperature  gradient  in  the  cell  as  the  result  of  the  increased  heat 
generation  rate.  However,  because  the  thickness  of  the  mem¬ 
brane,  CGDL  and  CCC  is  relatively  small,  the  overall  variation 
of  the  cell  temperature  in  the  cell  is  small. 

3.3.  Effect  of  the  channel  and  rib  width 

We  now  turn  our  attention  to  the  effect  of  the  channel  and  rib 
width.  The  simulations  were  performed  by  keeping  the  ratio  of 


DMFC. 
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Fig.  10.  Variations  in  the  operating  temperature  at  the  ACL  with  current  density 
for  different  channel  and  rib  width. 

the  channel  width  to  rib  width  at  1 : 1 .  Hence,  hereafter  a  specified 
rib  width  also  means  the  same  size  in  the  channel  width.  Also 
note  that  the  both  current  collectors  on  the  anode  and  cathode 
are  identical.  Figs.  9  and  10  show  the  cell  performance  and  vari¬ 
ation  in  the  cell  temperature  at  the  ACL  with  the  current  density 
at  2.0  M  methanol  operation,  respectively.  It  is  seen  from  Fig.  9 
that  the  cell  performance  decreases  with  the  rib  width.  This  is 
because  the  distance  of  mass  transfer  from  the  channel  region  to 
the  rib  region  increases  with  the  rib  width,  leading  to  a  higher 
mass  transfer  resistance  and  thereby  a  lower  methanol  and  oxy¬ 


gen  concentration  under  the  ribs.  Moreover,  the  larger  rib  width 
also  exacerbates  the  removal  of  the  gas  CO2  and  liquid  water, 
which  further  reduces  the  methanol  concentration  at  the  ACL 
and  oxygen  concentration  at  the  CCL.  As  a  result,  the  cell  per¬ 
formance  decreases  with  the  rib  width.  In  the  meantime,  at  a 
given  current  density,  the  overpotential  increases  with  the  rib 
width  as  the  result  of  the  increased  mass  transfer  resistance, 
resulting  in  a  higher  heat  generation  rate  and  thereby  a  higher 
cell  temperature,  as  shown  in  Fig.  10.  Particularly,  because  the 
effect  of  the  mass  transfer  becomes  more  significant  at  high  cur¬ 
rent  densities,  more  heat  due  to  the  mass  transfer  polarization 
is  generated.  Therefore,  the  temperature  difference  between  dif¬ 
ferent  rib  widths  increases  with  the  current  density.  However, 
although  the  cell  temperature  is  higher  with  the  larger  rib  width, 
the  cell  performance  is  still  lower.  This  fact  indicates  that  the 
effect  of  the  increased  mass  transfer  as  the  result  of  decreasing 
the  rib  width  is  more  important  than  the  effect  of  the  increased 
temperature  on  the  cell  performance.  Hence,  the  improved  cell 
performance  with  decreasing  the  rib  width  is  mainly  attributed 
to  the  increased  mass  transfer  rate  on  the  anode  and  cathode. 

Fig.  1 1  shows  the  distributions  of  methanol  concentration  at 
the  ACL,  current  density,  gas  saturation  at  the  ACL  and  liquid 
saturation  at  the  CCL  for  different  rib  widths  at  a  given  anode 
overpotential  of  0.5  V.  It  should  be  noted  that  the  rib  width  is 
normalized  to  the  total  width  of  the  channel  and  rib  for  the  pur¬ 
pose  of  comparison.  It  is  seen  that  the  smaller  rib  width  shows 
a  higher  methanol  concentration  and  more  uniform  distribution. 


Fig.  11.  Distribution  of  (a)  methanol  concentration  at  the  ACL,  (b)  current  density,  (c)  gas  saturation  at  the  ACL,  (d)  liquid  saturation  at  the  CCL  with  different 
channel  and  rib  width. 
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Fig.  12.  Distribution  of  (a)  cell  temperature  in  y  direction,  (b)  cell  temperature 
in  x  direction  with  different  channel  and  rib  width. 

This  is  because  the  smaller  rib  width  leads  to  a  shorter  distance 
of  mass  transfer  and  thereby  a  lower  overall  mass  transfer  resis¬ 
tance.  The  lowered  overall  mass  transfer  resistance  leads  to  a 
higher  methanol  and  oxygen  concentration  to  produce  a  higher 
and  more  uniform  current  density,  as  shown  in  Fig.  1  lb.  This  is 
also  the  reason  why  the  smaller  rib  width  yields  better  cell  per¬ 
formance  shown  in  Fig.  9.  As  mentioned  in  the  above  subsection, 
a  larger  rib  width  leads  to  a  larger  mass  transfer  resistance  of 
the  gas  CO2  from  the  rib  region  to  the  outside.  Hence,  the  gas 
saturation  under  the  rib  region  is  higher,  as  shown  in  Fig.  11c. 
Additionally,  although  the  generation  rate  of  the  gas  CO2  is 
larger  at  the  ACL  for  the  smaller  rib  width  because  of  the  larger 
current  density,  the  gas  saturation  at  the  ACL  is  still  lower  than 
that  of  the  larger  rib  width  as  the  result  of  the  lowered  mass 
transfer  resistance.  On  the  contrary  to  the  anode,  the  smaller  rib 
width  yields  higher  liquid  saturation  at  the  CCL  because  of  the 
higher  water  generation  rate  resulting  from  the  higher  current 
density,  as  shown  in  Fig.  lid.  Moreover,  the  liquid  saturation 
distribution  is  more  uniform  for  the  smaller  rib  width  as  a  result 
of  the  lowered  mass  transfer  resistance. 

The  temperature  distributions  in  the  v  and  y  directions  for 
different  rib  widths  at  a  given  anode  overpotential  of  0.5  V  are 
presented  in  Fig.  12.  From  this  figure,  it  can  be  seen  that  the 
larger  rib  width  exhibits  a  higher  cell  temperature.  This  can  be 
attributed  to  the  lower  water  evaporation  rate  resulting  from  the 
larger  rib  width.  Moreover,  the  current  density  at  the  channel 
region  is  higher,  resulting  in  a  higher  heat  generation  rate  and 
thereby  a  higher  cell  temperature  at  the  channel  region.  In  addi¬ 


tion,  the  cell  temperature  decreases  from  the  anode  to  the  outer 
surface  of  the  CCC  because  all  the  heat  is  released  to  the  ambient 
from  the  cathode.  It  is  interesting  to  notice  that  although  the  cell 
temperature  is  higher  with  the  larger  rib  width,  the  cell  perfor¬ 
mance  is  still  lower.  This  fact  reveals  that  the  mass  transfer  is 
a  predominant  factor  that  affects  the  cell  performance  when  the 
rib  width  is  varied. 

4.  Concluding  remarks 

A  two-dimensional,  two-phase  mass  transport  model  for  a 
passive  DMFC  with  the  heat  transfer  effect  has  been  developed 
based  on  the  unsaturated  flow  theory  in  porous  media.  A  com¬ 
puter  code  was  then  developed  for  solving  the  equations  that 
govern  the  interrelated  heat  and  mass  transport,  along  with  the 
electrochemical  reactions  occurring  in  the  passive  DMFC.  The 
effects  of  the  open  ratio  and  channel  and  rib  width  of  the  cur¬ 
rent  collectors  and  methanol  concentration  on  cell  performance 
were  investigated.  It  is  revealed  that  the  improved  performance 
as  the  result  of  higher  methanol  concentrations  is  due  primarily 
to  the  increased  cell  operating  temperature  resulting  from  the 
exothermic  reaction  between  permeated  methanol  and  oxygen 
on  the  cathode  and  increased  mass  transfer  rate  of  methanol.  It 
is  also  found  that  the  cell  performance  improves  with  increasing 
the  open  ratio,  which  results  in  the  improved  mass  transport  of 
both  methanol  and  oxygen,  although  the  cell  operating  tempera¬ 
ture  decreases  slightly  with  the  open  ratio.  Finally,  the  numerical 
results  show  that  the  cell  performance  degrades  with  increasing 
the  rib  width  (or  channel  width),  which  is  attributed  mainly  to 
the  increased  mass  transfer  resistance  at  the  both  anode  and  cath¬ 
ode,  although  the  cell  operating  temperature  increases  with  the 
rib  width. 
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